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Abstract. A crucial diagnostic of the ACDM cosmological model is the integrated Sachs- 
Wolfe (ISW) effect of large-scale structure on the cosmic microwave background (CMB). 
The ISW imprint of superstructures of size ~ 100 /i _1 Mpc at redshift z ~ 0.5 has been 
detected with > Act significance, however it has been noted that the signal is much larger 
than expected. We revisit the calculation using linear theory predictions in ACDM cosmology 
for the number density of superstructures and their radial density profile, and take possible 
selection effects into account. While our expected signal is larger than previous estimates, 
it is still inconsistent by > 3a with the observation. If the observed signal is indeed due to 
the ISW effect then huge, extremely underdense voids are far more common in the observed 
universe than predicted by ACDM. 
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1 Introduction 

The standard 'concordance' ACDM cosmological model fits many different observations, in- 
cluding the luminosity distance-redshift relation of Type la supernovae (SNe la) (e.g. [1-3]), 
the anisotropies in the CMB [4], the locally measured Hubble parameter (e.g. [5]) and baryon 
acoustic oscillations (e.g. [6] — but see [7, 8]). These observations, when interpreted assuming 
the homogeneous and isotropic Friedman- Robertson- Walker metric, imply that the expan- 
sion of the universe is accelerating, from which it is inferred that the universe is presently 
dominated by a cosmological constant ('dark energy') with negative pressure. It is important 
to note that this evidence is mostly geometrical, being based on interpreting measurements of 
distances — made using 'standard rulers' (the sound horizon at last scattering) and 'standard 
candles' (SNe la) — as due to accelerated expansion. The same data can equally be fitted 
without dark energy if, e.g., the isotropic but radially inhomogeneous Lemaitre-Tolman-Bondi 
metric is assumed and other assumptions such as a power-law spectrum for the primordial 
density perturbations are relaxed (e.g. [9, 10]). 

Given that dark energy is a complete mystery from a physical viewpoint, it is therefore 
imperative to examine the observational evidence for its dynamical effects. For example, the 
decay of gravitational potentials after dark energy begins to dominate (at redshift z < 1) 
should lead to secondary CMB anisotropies as the CMB photons traverse regions of over- or 
under-density — the ISW effect [11]. If the universe is spatially flat, then detection of the 
ISW effect through cross-correlation of the CMB with large-scale structure would provide 
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direct evidence of dark energy's negative pressure, hence crucial confirmation of the ACDM 
model [12]. 1 

To detect the ISW effect with 5a significance in CMB-galaxy cross-correlations requires 
z measurements for over 10 million galaxies [14, 15]. Such datasets are not yet available 
but several authors (e.g., [16-23]) have examined smaller source catalogues and reported 
detections with < 3a significance; however, [24] provides a skeptical view of some of these 
analyses. Conversely, some authors [25, 26] were unable to reject the null hypothesis (no 
ISW effect) and others [27] even found a slight anii-correlation rejecting ACDM at 2 — 
3a significance. Some groups have combined different data sets to increase the detection 
significance above 4<t [28, 29] but it has been argued that these analyses have underestimated 
the error bars [30]. 

Much of the uncertainty in full-sky studies arises from the difficulty in reconstructing 
the underlying density field from galaxy survey data, given Poisson noise in the galaxy 
distribution. A different approach to this problem is followed by Granett et al. [31, 32], 
who study the Sloan Digital Sky Survey (SDSS) Data Release 6 (DR6) luminous red galaxies 
(LRGs). They use 3D galaxy information rather than the projected 2D density, and select only 
the most extreme density perturbations, which are unambiguously identified despite Poisson 
noise. Along the lines of sight corresponding to these 'superstructures' they report a 4.4<r 
detection of the ISW effect, which is the most significant reported detection to date. This 
approach also provides information about the sizes and distribution of extreme structures 
in the universe so can be used to check the consistency of the standard ACDM model of 
structure formation, in particular whether the primordial fluctuations were indeed gaussian. 

However the magnitude of the temperature signal reported in Ref. [31] is surprisingly 
large and has been argued [33, 34] to be quite inconsistent with ACDM. In response, it 
has been noted [35] that the assumed profile of the superstructures has a big effect on the 
signal, and that the 'compensated top-hat' profile adopted in Refs. [33, 34] is not the most 
appropriate. Using an alternative profile that is calibrated against A^-body simulations and 
adopting a template-fit approach, it is claimed [36] that the discrepancy with ACDM is only 
at the 2a level. 

Our aim in this paper is to clarify this important issue. We calculate the expected 
temperature signal from these superstructures making no a priori assumptions about their 
nature except that they arose in a ACDM cosmology with gaussian primordial density pertur- 
bations. Our analysis uses linear perturbation theory, but non-linear corrections are argued 
to be sub-dominant at the relevant length scales and redshift. We use a statistical treatment 
of (initially) gaussian perturbations [37] to calculate the expected number densities and the 
expected profiles of the superstructures. These profiles are not significantly different from 
those used in Ref. [35] , and we do find a factor of ~ 3 increase in the ISW signal over that cal- 
culated earlier [33, 34]. However, contrary to Ref. [36] we find that the expected value of the 
signal is still discrepant at > 3a with the observations reported by G08a. We demonstrate 
that this difference in the conclusions arises due to the interpretation of the template-fitting 
approach used in Ref. [36] and argue that a correct interpretation would lead to conclusions 
compatible with those presented here. 

In Section 2 we briefly review the ISW effect and in Section 3 we calculate the expected 
temperature signal of superstructures in the standard ACDM model. In Section 4.1 we 

1 The ISW effect should also boost low multipoles in the CMB angular power spectrum, whereas these are 
in fact anomalously low on the observed sky. However given the large cosmic variance on these scales and 
galactic foreground systematics, this discrepancy is not thought to be significant [13]. 
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describe the key features of the observational strategy [31, 32] which must be accounted for 
before making a comparison with the theoretical calculation. In Section 4.2 we show that 
even if Refs. [31, 32] had selectively picked out the regions in the survey with the biggest ISW 
signal, there is still a significant discrepancy with the ACDM expectation. In Section 4.3 we 
discuss the reasons for the difference between our result and that of Ref . [36] . 

Finally in Section 5 we turn to possible alternative explanations outside the standard 
cosmological model. Modifications to the growth rate of perturbations from the ACDM 
prediction could be responsible for the discrepancy, but any alternative model must also 
match the other observational successes of ACDM. A plausible explanation is that the effect 
is due to non-gaussianity of the primordial perturbations, but further work is needed to 
calculate its effect on the expected signal. Until this discrepancy is clarified, the observations 
reported in Ref. [31] do not appear to be compatible with the ISW effect expected in the 
standard cosmology. 



2 The ISW effect 



In a universe with matter density £l m = 1 and no dark energy, density perturbations 5 
grow at exactly the same rate as the scale factor of the universe a (= (1 + z)~ l ), so at 
the linear level there is no evolution of the gravitational potential (oc —8/a). However, 
in a ACDM universe, a grows faster than (linear) density perturbations, so perturbations 
in $ decay with time. For a CMB photon passing through an overdense region the energy 
gained while falling in is not cancelled by the energy lost in climbing out of the evolved, 
shallower, potential well. Overdense regions (clusters) therefore appear as hot spots in the 
CMB; conversely, underdense regions (voids) will appear as cold spots as the photon loses 
more energy climbing the potential hill than it gains subsequently while descending. 
The temperature fluctuation AT(n) induced along direction n is [11]: 

2 - r L ■ 

AT(n) = -g-Tb / ®(r,z,n)adr, (2.1) 
c Jo 

where To is the mean CMB temperature, tl is the radial comoving distance to the last 
scattering surface (LSS), $(r, z, h) is the time derivative of the gravitational potential along 
the photon geodesic and c is the speed of light. 

The Poisson equation relates <I> to the density contrast 5 = (p — p)/p (where p(t) is the 
mean density) through: 

V 2 $(x, t) = 4vrGp(t)a 2 (5(x, t). (2.2) 
This can be written in Fourier space as 



where Hq is the current Hubble parameter. Taking the time derivative of this equation yields: 

(2.4) 



in 



a 



We assume that linear theory holds on the large scales of interest hence perturbations 
grow as S(k, t) = D(t)S(k,z = 0), where D(t) is the linear growth factor. A numerical 
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simulation has shown that non-linear effects represent only a 10% correction at the low 
redshifts we are interested [38]. 2 In this approximation, 
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where f3(z) = dlnD/dlna is the linear growth rate. Hence the time evolution is captured 
by the ISW linear growth factor, G(z) = H(z) (1 — f3(z)) D(z)/a. For an J7 m = 1 universe, 
P(z) = 1 for all z so there is no ISW effect. 

Given the density profile b of any isolated superstructure, Eqs. (2.1) and (2.5) can be 
used to calculate the temperature fluctuation it induces in the CMB. Assuming spherical 
symmetry of the density profile, Eq. (2.5) in real-space becomes: 



with S(r') evaluated at redshift z = 0. Thus F(r) contains all information about the structure 
in question, while the assumed cosmology enters through the prefactor and the ISW growth 
factor G{z) in Eq. (2.6). 

3 Expected signal from superstructures in ACDM 

'Superstructures' refer to density perturbations extending over > 100 h -1 Mpc and should 
not be thought of as non-linear collapsed structures in the usual sense, rather as smooth hills 
and valleys in the density distribution. Collapsed structures form only where the density 
perturbation 5(r) exceeds unity, which happens on much smaller scales than those of interest 
here. 

It is stated in Ref. [31] that the most extreme structure in the (500 /i _1 Mpc) 3 box of 
the Millennium A-body simulation [39], when placed at z = 0, would produce a signal of 
AT ~ 4.2 /xK. However when the signal distribution for supervoids with the densities and 
sizes reported in Ref. [32] was calculated, the answer was only (AT) = —0.42 /iK [33]. These 
authors [33] assumed a 'compensated top-hat' density profile motivated by the asymptotic 
final state of a void [40]. A similar profile was assumed by other authors [34] who found 
a similar average signal (AT) = —0.51 /uK for the 50 most extreme density perturbations 
of fixed radius r = 130 /i _1 Mpc expected in a ACDM cosmology. Subsequently it was 
argued [35, 36] that this profile is not the appropriate choice for density perturbations on 
> 100 /i _1 Mpc scales; instead an uncompensated density profile was motivated from gaussian 
statistics, and was found [35] to give larger values of AT. Using this profile Ref. [36] claimed 
only a 2a discrepancy between the ACDM prediction and observation. However, as we discuss 
in Section 4.3, the interpretation of [36] implies underdense regions with a physical density 
contrast 5 < — 1, which is physically impossible. 

It is thus necessary to revisit this issue. Using the statistics of a homogeneous, isotropic, 
gaussian density field, we now derive the expected mean density profiles of superstructures 

2 Both the linear and non-linear effects grow with time, however at late times and large scales (~ 
100 /i _1 Mpc) the linear effect dominates while at early times (when Q.a — 0), both effects are smaller but the 
non-linear effect dominates. 
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of all density contrasts and all sizes, as well as the expected number density of such super- 
structures. 



3.1 The number density of structures on different scales 

We identify superstructures of different sizes with extrema of the linear density perturbation 
field 5(r) when smoothed over different scales. Overdensities correspond to peaks of the 
smoothed field and underdensities to troughs. In the ACDM model, 5(r) is a homogeneous 
and isotropic, gaussian-distributed random field and the statistical properties of the maxima 
and minima have been calculated in Ref. [37]. We briefly review below their key results and 
introduce necessary notation. 

Let P(k, t) denote the matter power spectrum, defined as the Fourier transform of the 
2-point correlation function £(r, t) of the density field at time t. Define a set of spectral 
moments weighted by powers of k: 

/h 2 c\k 
—W 2 (kRf)P(k,t)k 2 \ (3.1) 

where W(kR{) is the window function appropriate to the filter used to smooth the density 
field, and -Rf is the (comoving) smoothing scale. Thus <jq is just the standard rms fluctuation 
of the smoothed density field. Using a gaussian filter, W(kR{) = ex.p(—k 2 R 2 /2), we define 
the parameters: 

7=^, R. = y/g±. (3.2) 

0-200 £T 2 

The (comoving) differential number density A/ r max ( z/ ) of maxima of height 5q = va§ is then 
[37]: 

M^{v)dv = (2 ^| e-^/ 2 G( 7 , iv)&v • (3.3) 

The function (7(7,71/) is given by Eq. (A19) of Ref. [37]; we use a fitting form, accurate to 
better than 1%, given in their Eqs. (4.4) and (4.5). The density of minima is related simply 
to that of maxima through: N" m i n (f) = Nmax(— l/ ). 

3.2 Mean radial profiles 

Having identified superstructures with the maxima or minima in the smoothed density field, 
we wish to determine the mean radial variation of the density field in the neighbourhood of 
these extrema. Ref. [37] shows that, given a maximum 5 = 5$ at r = 0, the mean shape in 
the vicinity of this point after averaging over all possible orientations of the principal axes 
as well as all values of the curvature at r = is: 

*M = (* + ? vv )- (tV + • (3-4) 



(1-7 2 ) V 3 / 7(1"7 2 ) V 3 
where ifj(r) = £( r )/£(0) is the normalised density-density correlation function, and x = 
— V 2 <5/o"2. The expectation value of x given a peak of height 5q at r = is approximately 

(x\S ) =71/ + 9(7,7!/) , (3-5) 

where v = Sq/gq, and 0(7,71/) is given by the fitting function (6.14) of Ref. [37]. It follows 
that Eq. (3.4) can be rewritten [41] as: 



5(r,t) = if ^^f^W\k Ri )P(k,t) 
do Jo 2tt z kr 
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dJfe, (3.6) 
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Figure 1. Left panel: Mean radial profiles for voids obtained from Eq. (3.6) for different values of 
central underdensity So and the smoothing scale R[. The blue (solid) line and the green (dash-dot) 
line are for the same Rf (=20 ft. _1 Mpc) but different values of So, whereas the red (dotted) profile 
has the same <5o (= —0.2) as the blue (solid) profiles, but a larger smoothing scale. The blue (dashed) 
curve is the biased galaxy density contrast corresponding to matter density contrast given by the blue 
(solid) line, with bias factor b = 2.25 as is appropriate for LRGs. Right panel: Comparison between 
the two methods of calculating mean radial profiles. The black (dash-dot) line represents <5 r andom(7") 
from Eq. (3.7) for So = —0.3 and Rf = 30 h~ 1 Mpc. The blue (solid) line is the corrected profile S(r) 
from Eq. (3.6) for the same parameter values. 



which we use for numerical evaluation of profiles. Note that 7 and i?* depend on the smooth- 
ing scale Rf. 

We make the simplifying approximation that the average ISW signal for a large number 
of superstructures is the same as the signal due to a superstructure with the mean profile 
(3.6) — this is justified because of the linear relationship between AT, and 5. The actual 
distribution of the size of structures depends on the distribution of x values, which can be 
obtained from Eq. (7.5) of Ref. [37]. In the left panel of Fig. 1 we plot some underdense 
profiles for selected values of 5q and -Rf (these examples are chosen for clarity and are not 
representative of the most likely actual underdensities). Note that the size of the structures 
is much larger than the smoothing scale -Rf. 

To identify superstructures in galaxy surveys (the methodology of Ref. [31] is discussed 
in more detail in Section 4.2) the galaxy density contrast 5 g is assumed to be linearly biased 
with respect to the matter density: 5 g = b5. Denoting by p s \ the value of the density 
field at turnover, and by po the minimum density at the centre, a selection cut is made 
on on w = /9 s i/pcb amounting to a lower bound on the absolute value of 5q. This avoids 
false detections of over- and under-dense regions that are just Poisson fluctuations in galaxy 
number counts and ensures that only extreme superstructures are included in the ensemble. 

The definition of the radius R v of any of the voids shown in Fig. 1 is slightly ambiguous. 
We choose it to be the radius of turnover in the density profile less Rf, since smoothing 
necessarily increases the radius somewhat. This is a small correction since in general R v S> 
R f . 

The mean profile about a random point with a density contrast 5q (i.e., a point which 
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is not necessarily an extremum of the density field) is simply [42] 



^random(^) = 5^{r) . (3.7) 

This should be compared with Eq. (3.4), which is the appropriate mean profile about points 
that are also known to be extrema. Eq. (3.7) is the form of the density profile used in 
Refs. [35, 36]. In Ref. [36] this is matched to the mean density profile about randomly chosen 
points in TV-body simulations and is found to be a good approximation, as expected. It must 
be emphasised that as we identify superstructures with points of extrema in the density 
field, consistency demands that we use the corrected form of Eqs. (3.4) and (3.6) rather than 
Eq. (3.7). This gives a profile that is slightly steeper than that used in [35, 36]. 

In practical terms, however, the choice of profile between Eq. (3.4) and Eq. (3.7) will 
have little effect on our conclusions as the difference between these profiles is small. In the 
right panel of Fig. 1 we show as an example the mean profiles 6(r) and <5 ran dom for 5o = —0.3 
and R t = mh^M^c. It can be seen that the two are similar: in fact, for most values of 
5q and Rf and over almost the entire range of r, 5 r andom( r ) lies within the ±lcr dispersion 
around 5(r) of the profiles of extrema (this is not shown here but some examples can be seen 
in Fig. 8 of Ref. [37]). 

Of course the quantity determining AT in Eq. (2.1) is not 5(r) but an integral over F(r), 
which is itself an integral over 5(r), and hence small differences may be magnified. Therefore, 
to explicitly check that the choice of profile is not responsible for the difference between our 
results and those of [36], we also repeat the calculation of the expected temperature signal 
described below using <5 r andom(?") to model structures. 

3.3 Temperature signal 

The ISW signal of any individual superstructure will be too small compared to the primordial 
CMB anisotropies to be observable. Therefore, what is measured is the average temperature 
fluctuation along the lines of sight of a selected sample of either over- or under-densities. The 
primordial anisotropies are uncorrelated with the large scale structure and average out so 
given a large enough sample the correlated ISW signal eventually dominates. Our calculation 
of this averaged signal is done as follows. 

At a given z, we use Eq. (3.6) to calculate the matter and galaxy density profiles about 
extrema of the density field as functions of 5q and Rf and obtain along the line of sight 
as discussed in Section 2. This enables us to calculate AT(8; 5q, Rf) where 9 = 0° is the line 
of sight passing through the centre of the superstructure. To compare with the observations 
we first apply the selection criterion on 5q through the limit on w. Then, to calculate the 
expectation value (AT) for the resulting ensemble, we weight the results appropriately with 
the number density (3.3) of extrema. Hence for an ensemble of voids: 



// W( e)AT(9-5 )SS min cTo l d 2 6dS 
7r6>2 / M min a 1 d5 



(AT) = JJ " ^^^^"o - ^" u , (3i 



where < 6 < O ut] W(0) is a filter of angular width 9 ou t chosen in order to match that 
used in the actual observation and —1 < 5o < 5q, where <5§ is the (mildly i?f-dependent) 
cutoff on the minimum underdensity required to pass the selection criterion. The choice of 
i?f determines the mean radial size of the voids included in the ensemble; although structure- 
finding algorithms may not have an explicit size dependence, in practice there is obviously a 
lower limit on the size of the over- or under-density that can be reliably found. As smaller 
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structures are overwhelmingly more probable than larger ones, it is important to capture 
this effect and we discuss this in more detail in the next Section. For a given ensemble of 
voids specified by Rf and Sq, the observed temperature signals AT will be distributed about 
this mean value. By calculating (AT 2 ) in a similar fashion to Eq. (3.8) we can estimate the 

1 /2 

standard deviation of this distribution simply as ((AT 2 ) — (AT) 2 ) 

The expected signal from an ensemble of clusters follows in an exactly analogous manner 
to Eq. (3.8). 

4 Comparing theory to observation 

4.1 The measured ISW signal of superstructures 

To compare our expectation for the ISW signal to the measurement made by Ref. [31], it 
is necessary to follow the same selection procedure. They use a sample of 1.1 million LRGs 
in the range 0.4 < z < 0.75 (median z = 0.52) from the SDSS DR6 [43], which covers 
7500 degree 2 on the sky and occupies a volume of 5 /i~ 3 Gpc 3 . They search for 'supervoids' 
and 'superclusters' using two publicly- available structure-finding algorithms: ZOBOV (ZOnes 
Bordering On Voidness; [44]) for supervoids, and VOBOZ (VOronoi BOund Zones; [45]) for 
superclusters. 

It is necessary to mimic the way in which these algorithms select structures in choosing 
the ensemble for which to calculate (AT) from Eq. (3.8). ZOBOV uses a parameter-free 
Voronoi tessellation to estimate the density at each galaxy in the sample, based on the 
distance to its nearest neighbours. Around each density minimum it then finds the region 
of the density depression or supervoid. (Of course large voids can contain multiple smaller 
voids, or even isolated high-density regions.) The 'significance' of the depression is estimated 
by comparing the density contrast, w (defined as the ratio of the density at the lip of the 
void to the density at its minimum) to an uniform Poisson point sample. This yields the 
likelihood that a void of density contrast w could arise from Poisson noise, i.e. that it is 
a false positive detection; a 3a cut is then applied on the likelihood which translates to 
requiring w > w c = 2.0 on the density contrast [44]. This procedure yields 50 supervoids, the 
properties of which are tabulated in Ref. [32]. The VOBOZ supercluster finder uses the same 
algorithm but applied to the inverse of the density field, with density contrast defined as the 
ratio of the peak density to the density at the edge of the structure. However, now the 3a 
cut on the likelihood that an overdensity of given w could have arisen due to Poisson noise 
corresponds to w > 6.8 [45]. In fact the authors impose a tighter cut: w > 8.35 in order to 
obtain exactly 50 such superclusters; their properties are also tabulated in Ref. [32]. 

Ref. [31] report a search for the ISW signals of these superstructures using an inverse- 
variance weighted combination of the WMAP 5-year Q, V and W maps [46] with foreground 
subtracted and the KQ75 mask applied. They build stacked images by averaging the CMB 
temperature in the regions around the lines of sight passing through the centres of the 
identified superstructures and use a compensated top-hat filter of width 9 C in order to perform 
the averaging. This corresponds to making the choice in Eq. (3.8): 

W{9) = { lj °^ e ^ 9 » (4.1) 
[ —1 , 9 C < 9 < #out , 

with 6> ou t = V29 c . For a filter radius 9 C = 4°, 3 the sample of supervoids gives (AT) = 
— 11.3 ± 3.1 //K, and the sample of superclusters (AT) = 7.9 ±3.1 fiK. When averaged 

3 This choice is apparently motivated by the expectation that the CMB-galaxy cross-correlation should 



- 8 - 



together (with the negative of the supervoid image added to the superclusters) this gives 
(AT) = —9.6 ± 2.2 /jK, i.e. a 4.4<r detection of a non-zero correlation. 

The quoted error bars are determined by Monte Carlo simulations in two different 
ways: by stacking images of the actual CMB map in random directions, and by keeping the 
directions fixed but using model CMB maps. Note that this explicitly accounts for cosmic 
variance in the CMB fluctuations. The theoretical error bars described in Section 3.3 account 
for the distribution of voids of different depths and sizes. Misidentification of structures in 
the dark matter density field can also occur due to shot noise in the galaxy distribution. 
This means that the structures identified as the N most extreme fluctuations will in fact be 
a random sample of N structures drawn from among the M most extreme, for some M > N. 
In Section 4.2.1, we make the optimistic assumption that N = M such that the expected 
signal is maximised. 

A caveat that should be noted, however, is that the choice in Ref. [31] of using exactly 
50 supervoids and superclusters each is made to obtain the largest signal-to-noise ratio. For 
both A = 30 and N = 70 a somewhat lower significance detection is reported compared to 
N = 50. As the distribution of AT values for the voids is very non-Gaussian, this choice of ./V 
could be important. Fluctuations due to the underlying CMB anisotropy dominate at small 
A and false positive identifications of structures increase at large A so the signal-to-noise 
ratio is indeed expected to have a maximum at some intermediate A. However more detailed 
study is needed to quantify whether this is expected to be at A ~ 50. In this paper we 
take the reported observations at face value and our interpretation is based on the data for 
A = 50. 

4.2 Comparison to the theoretical expectation 

We assume for simplicity that all the superstructures are located at the mean redshift z = 0.52 
and adopt the standard ACDM cosmological model with Q m = 0.29, = 0.71, n s = 0.96 
and eg = 0.83 (with h = 0.69 where required) to obtain the matter power spectrum at 
z = 0.52 using CAMB [47]. 4 The bias factor for LRGs is taken to be b = 2.25. These 
parameters are kept fixed for the purposes of calculation, but as we discuss later, varying 
them within the bounds set by other observational data has no effect on the conclusions. 

Our first finding is that there are no overdense superstructures within the linear regime 
(i.e. with 5q < 1) which meet the VOBOZ 3<7-significance selection criterion that w > 6.82. 
Such a ratio of densities (between the lip of an overdensity and its centre) can be achieved, 
but only for non-linear collapsed structures. We conclude that Table 5 of Ref. [32] does 
not list the most overdense large-scale linear perturbations, but the (mild) large-scale linear 
perturbations that happen to contain the most overdense small-scale collapsed structures. 
This means that the relative ranking and selection of large-scale linear overdensities made 
in Ref. [32] is done on the basis of the small-scale collapsed structures they contain. The 
statistical effect on the CMB temperature of small collapsed structures could be calculated, 
but we cannot model the criteria by which these overdense regions were selected for study. 
Note however that our calculations below for the maximum possible amplitude of ISW signal 
from superstructures holds equally well for over- and under-densities in the linear regime. 

peak at about 4° [22]. Ref. [31] reports repeating the observation with a few other widths, 3° < 6 C < 5° and 
obtaining a maximum detection significance for C = 4°. 

4 These are the mean parameter values obtained from a fit to WMAP 7-year [4] and SDSS DR7 [48] data 
using COSMOMC [49]. 
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Any contamination of the regions selected by VOBOZ will only reduce the expected signal, 
although we are unable to estimate by how much. 

On the other hand, this problem of modelling the selection of structures for study does 
not arise for voids. Hence we concentrate only on the sample of underdense regions (i.e. 
supervoids), for which (AT) b s = —11.3 ± 3.1 fiK. To calculate the expected (AT) using 
Eq. (3.8), we must first choose the smoothing scale Rf and determine the distribution of void 
radii R v in the ensemble using Eq. (7.5) of Ref. [37]. For each Rf we calculate i?™ m such 
that 95% of all voids in the ensemble have radius R v > R^ nn , this being a convenient way 
to characterise the ensemble. For Rf = 20 /i _1 Mpc, we find i?™ in ~ 70 h M.pc, which is 
similar to the mean radius of the voids in Table 4 of Ref. [32] ; 5 at this scale, the number 
of such voids that should satisfy the selection criterion is N v ~ 10 4 . This is to be compared 
with the N v = 50 voids that are actually tabulated in Ref. [32] . For the larger ensemble, we 
find an expectation value (AT) = —0.3 ± 0.2 fiK, i.e. consistent with zero. 

However, since Granett et al. see only a tiny fraction of the total number of supervoids 
expected in the SDSS DR6 volume, it is unreasonable to assume that this represents a fair 
sample. Some strong selection effects must be in operation; while these have no effect on our 
ISW prediction for individual voids, they will affect the selection of structures in the LRG 
distribution and thus the ensemble expectation value. Selection effects that can enhance the 
expected signal are a bias towards larger and deeper regions, implying that Granett et al. 
did not randomly select 50 of the ~ 10 4 expected supervoids but chose some sample that 
is skewed towards regions with larger AT values. We show below that the expected signal 
from the 50 most extreme regions is indeed ~ 5 times larger than —0.3 [aK; nevertheless the 
discrepancy with the observation in Ref. [31] is still > 3a. Therefore, irrespective of how the 
regions were selected, tension remains with the expectation in the standard ACDM model. 

4.2.1 Accounting for selection effects 

We first consider the possibility that the void-finding algorithm ZOBOV is sensitive only to 
the largest (and least common) voids in the matter distribution which produce the biggest 
ISW temperature signals. The average density of LRGs in the SDSS DR6 sample is roughly 1 
galaxy per (15 /i _1 Mpc) 3 . In underdense regions LRGs will be even more sparsely distributed 
so ZOBOV will certainly be less able to identify smaller underdense structures, thus biasing 
the sample towards larger voids. 

In order to model the effect of such a selection bias, we increase the value of Rf in 
Eq. (3.8); this is equivalent to including only the N v largest voids with radius R v > i?™ m 
in the ensemble from which (AT) is calculated. In the left panel of Fig. 2 we plot the 
expectation (AT) as a function of R™ m . The orange cross-hatched area shows the region 
that is within 3a of the observed value (AT) t, s = —11.3 ± 3.1 /iK. The theoretical value 
of (AT) becomes marginally consistent with the observed value when i?™ m ~ 170 /i _1 Mpc. 
However the probability that within the SDSS survey volume there are 50 supervoids of 
radius R v > 170 /i -1 Mpc which also meet the ZOBOV selection criterion is vanishingly small. 

It is interesting to note that the 50 largest supervoids expected within the SDSS survey 
volume have i?™ m > 120 h M.pc. The largest void radius reported in Ref. [32] is 125 ft. Mpc 
and the mean is 70 /i _1 Mpc. We have argued that these values somewhat underestimate 
the size of the supervoids compared to our definition of R v , yet it seems unlikely that the 
difference could be so large that all the tabulated voids should have R v > 120 ft Mpc. 

5 At the radii reported here, the density has not yet reached the background level so the quoted values 
must be underestimates of the void size relative to our criterion. 
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Figure 2. Left panel: The absolute value of (AT) for an ensemble of voids which satisfy the 
ZOBOV selection condition on density (see text), as a function of the minimum radius of voids in the 
ensemble. The solid (blue) curve shows the mean value and the shaded (lighter blue) contours the 
la region. The (orange) cross-hatched area is the lower end of the 3er range of the observed value 
(Ar) bs = —11.3 ± 3.1 /iK. Right panel: As above, but showing (AT) as a function of the number of 
voids in the ensemble from which the observed sample of 50 voids is to be drawn, when only the N v 
largest voids also meeting the ZOBOV selection condition on density are included. 

In the right-hand panel of Fig. 2 we plot (AT) as a function of the size N v of the ensemble 
of the largest supervoids that should exist within the SDSS volume. It is from this ensemble 
that the 50 observed supervoids should be regarded as having been drawn. It can be seen 
that even under the assumption that the VOBOZ algorithm selected exactly the 50 largest 
voids in the entire SDSS survey volume, the expected signal is only (AT) = —1.33 ±0.13 /uK 
which is still discrepant by > 3<r with the observed value. We conclude that the observed 
signal cannot be explained due to a simple bias towards selecting only the largest voids. 

We repeated the calculation above using Eq. (3.7) to model the density profiles of the 
voids, as in Ref. [36], and obtained (AT) = —1.58 ±0.12 /iK for the 50 largest voids. Despite 
the small increase compared to the value quoted above, this is still discrepant with observation 
at the same > 3a level. This conclusively demonstrates that, as anticipated, the difference 
between our conclusion and that of Ref. [36] is not due to the small correction included in 
S(r). Henceforth we use the corrected profile only. 

We consider next whether the ZOBOV algorithm is more sensitive to deeper voids. In 
Table 4 of Ref. [32], the edge of most of the supervoids is defined at a radius where the 
density contrast is still negative. This means ZOBOV systematically underestimates the 
value of w relative to our definition (where S e d ge — 0), so Granett et al. effectively used a 
more stringent cut on w than we have done. We can model this effect by varying 5q from 
the value determined by the stated algorithm. In the left-hand panel of Fig. 3 we plot as 
examples (AT) as a function of 5q for i?™ m ~ 70 /i _1 Mpc (the mean radius of the supervoids 
in [31]) and i?™ m ~ 100 h~ 1 Mpc. The right-hand panel shows (AT) as a function of N v , 
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Figure 3. Left panel: The absolute value of (AT) for an ensemble of supervoids which satisfy Sq < Sq, 
as a function of Sq. The dashed (green) curve shows the case when Rf is chosen such that voids with 
R v > 70 /i _1 Mpc are included in the ensemble; the solid (blue) curve is for R v > 100 h~ 1 Mpc. Shaded 
contours show the la region about the mean and the orange cross-hatched area is as in Fig. 2. Middle 
panel: As before, but showing (AT) as a function of the number of voids in the ensemble from which 
the observed sample of 50 supervoids is to be drawn, when only the N v deepest supervoids are included. 
Right panel: As before, but for two different choices of the radius 9 C of the compensating top- hat filter 
used in Eq. (3.8). The solid (blue) curve is the mean value for 8 C =4°, and the dot-dashed (black) 
curve for 9 C = 6°. Shaded regions show the la deviations from the mean. 



the number of voids included in the ensemble when Sq is varied in each case. For the smaller 
radius, the ensemble is dominated by small voids with a small ISW effect, so increasing i?™ m 
increases (AT) at any Sq. With R™ m ~ 100 /i -1 Mpc, (AT) becomes marginally consistent 
with the observation for Sq < —0.5. However, the probability of obtaining 50 supervoids 
with R v > 100 /i^Mpc and S < -0.5 is negligibly small, as the right-hand panel clearly 
demonstrates. Therefore we conclude that in the standard cosmological picture there do not 
exist 50 isolated matter density voids in the SDSS volume that, irrespective of the method of 
selecting them, produce an average (AT) compatible with observation. Our analysis, based 
on very conservative assumptions, shows that the claim [36] that the discrepancy between 
observation and the ACDM prediction is only ~ 2a is an underestimate. 

4.2.2 Potential systematic errors 

We now discuss the expected corrections due to the simplifying assumptions we made in 
calculating the signal. We first consider the effect of varying the cosmological parameters 
from the mean values adopted in 4.2. This can affect the calculation of (AT) in three ways: 6 
by altering the mean radial profile of the voids; by altering the expected number density 
of voids through Eq. (3.3); and by altering the ISW growth factor G(z) and the integral 

6 Varying the bias factor only affects the detection of matter density voids via the observation of galaxies. 
By varying Sq as above we have already made the conservative assumption that only the most extreme matter 
voids that exist are detected. Therefore the bias factor can have no further effect on our conclusions. 
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Eq. (2.1). The first two of these are dependent, through the moments <Tj, on the matter power 
spectrum P{k). An increase in power increases the number density of points of extrema, and 
steepens the mean profile about these points. However, it also increases as, which is already 
reasonably well constrained by WMAP [4] and SDSS DR7 [48] data. Within any reasonable 
range of as, the change to the mean profile 5(r) is completely negligible. Allowing as = 0.87, 
which is the 95% C.L. bound from WMAP + SDSS, does boost the relative number density 
of the largest voids by a factor of ~ 1.5. However, given the almost logarithmic dependence 
of (AT) on N^ g evident from Fig. 2, the effect on (AT) remains well within the statistical 
error bars already quoted. 

The only remaining way for the cosmological parameters to affect (AT) is through the 
value of Q m . Assuming a spatially flat model, this determines the growth factor G(z), the 
constant pre-factor in $ and the photon geodesies (i.e. the relationship between r and z in 
Eq. (2.1)). Note that a smaller f2 m boosts G(z) but decreases the pre-factor. f2 m is of course 
also constrained by other data, and we find that any value in the range Q m £ (0.25, 0.32) 
results in changes to (AT) that are significantly smaller than the error bars already quoted. 
We conclude that the freedom to vary the ACDM model parameters is not enough to overcome 
the discrepancy with observation. 

We have neglected non-linear effects and the time evolution of the void density profile 
[40, 50] in the calculation of (AT). On small scales voids evolve non-linearly towards a 
compensated top-hat profile. These non-linear effects increase the depth of the cold spot 
at its centre but also produce a hot filamentary shell [38]. When averaged over angle 9 it 
is not clear whether the net effect will lead to a net increase or decrease in (AT). In fact, 
highly non-linear void profiles are known to produce a smaller ISW signal due to the effect 
of the overdense ridge at the boundary [34, 35], so by using the linear theory profile we may 
be overestimating the expected signal. It must be stressed though that at small redshifts 
z ~ 0.5 and on the large scales of interest > 100 /i _1 Mpc, the effects of non-linear evolution 
will in any case represent only a subdominant correction [38] and can justifiably be neglected 
here. 

A more subtle assumption is that the real supervoid profiles are adequately captured 
by our smoothing prescription. The effect of the smoothing is to slightly broaden the <5(r) 
profile which in turn leads to broadening of the AT (6) profile. This may become a problem 
if the AT (6) profile is significantly broadened so the compensated top-hat filter of radius 
4° in Eq. (3.8) then underestimates the real signal. If this were the case then it would be 
more appropriate to use a broader top-hat. To check this, we repeated our analysis with 
a filter of radius 6°, which is a generous overestimate of the degree of broadening, and in 
the right panel of Fig. 3 we show the effect for the largest and deepest supervoids. As 
expected, increasing the filter radius does increase (AT), but the effect is small. Even with 
a 6° filter the expected signal remains > 3a discrepant with observations for A^ v = 50; also 
as mentioned above, the actual effect on (AT) due to the smoothing will be less than this 
extreme model. This indicates that our conclusion of a discrepancy will be independent of 
the details of the smoothing filter used. 

We have made the simplifying assumption that the centres of all the superstructures 
lie at the same redshift z = 0.52, which is the median redshift of the DR6 survey. This 
means we have not accounted for the possibility that two voids may lie at different redshift 
along roughly the same line of sight, thus increasing the temperature shift. Equally we do 
not account for the possibility that a void may lie in front of or behind an overdensity, thus 
decreasing the temperature shift. We expect that these two simplifications largely negate 
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each other. 

As a test of the robustness of our calculation, we compare our results with Fig. 1 of 
Ref. [38] which shows the ISW map from a cosmological iV-body simulation for a volume 
comparable to the SDSS DR6, but at z = rather than z = 0.52. The very largest density 
perturbations in this map yields a maximum ISW signal of |AT| ~ 4 /iK before applying 
a filter analysis. Taking into account that AT is more pronounced at smaller redshift, this 
tallies very well with our prediction from Fig. 3 that the most extreme supervoid in the SDSS 
volume should produce (AT) ~ — 2 /xK. 

As seen in the full-sky maps in Ref. [38] , there are lines of sight along which the chance 
alignment of several large structures between us and the LSS several Gpc away can lead to 
isolated ISW 'cold spots' with AT < — 10 fiK. However, this is a cumulative effect of several 
structures over Gpc scales and does not affect our conclusions about the correlation between 
these cold spots and isolated structures at z ~ 0.5. (Note also that in the simulation the 
temperature profiles along such lines of sight are much broader that 4° and so when averaged 
using a 4° compensated top-hat filter would still give only a small |(AT)|.) Similarly, there 
will be several 4° circles on the sky for which the underlying CMB anisotropy alone can 
give an average AT of similar magnitude — it is just such fluctuations that generate the 
observational uncertainty of ±3.1 fiK of Ref. [31]. Such fluctuations are not due to any ISW 
effect and there is no reason why they should be correlated with large structures in the SDSS 
catalogue. Therefore, assuming only that the lines of sight were not chosen a posteriori, 
the observed correlation between galaxy structures and the CMB can only result from rare 
(> 3a) fluctuations or anomalously large density perturbations at z ~ 0.5. 

Ref. [31] also reports the variation of the signal-to-noise ratio with the width of the 
compensating top-hat filter and finds a maximum at 8 C = 4° for the combined sample of 
over-and under-dense structures. By contrast, for supervoids we find that |(AT)| increases 
as 8 C is raised from 4° to 6° (Fig. 3) and to larger angles. This is incompatible with a peak 
in the signal-to-noise ratio at 4° and is another respect in which the observation of Ref. [31] 
cannot be modelled by linear structures in ACDM. 

4.3 Comparison with earlier results 

We have provided a theoretical calculation of the expected average temperature signal in 
the CMB from superstructures observed in the SDSS LRG catalogue assuming a ACDM 
cosmology with gaussian primordial density perturbations. On comparing this expectation 
to the observation made in Ref. [31], we find a discrepancy of > 3a, even if possible unknown 
selection effects are accounted for. In fact we find that the single most extreme void expected 
within the SDSS survey volume does not produce an angle-averaged AT value within 3a of 
the observed average signal for 50 voids in G08a. On the other hand, in Ref. [36] Papai et al. 
perform a similar comparison using a different method, and report a discrepancy with ACDM 
of only ~ 2a. As demonstrated in Section 4.2.1, this is not due to the slightly different form 
of the radial profile for superstructures used. Instead the source of the difference arises from 
the interpretation of the matched filter technique employed in Ref. [36] . 

Papai et al. calibrate their radial profiles against simulation to obtain a template for 
the ISW effect from structures. As they note, reconstructing accurate estimates of the ISW 
profiles from comparison with simulations of finite size is complicated by cosmic variance in 
the simulation, necessitating the use of a high-pass filter in the reconstruction. Once this 
filter is applied, their numerical template profiles agree well with our analytic ones. 
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They then assume exactly the same template profile for all structures, scaled by choosing 
60 = 1 for superclusters and 60 = —1 for supervoids, and place these structures at redshift 
z = 0.52 at the locations of the structures observed in G08a. However, this procedure does 
not account for the actual density contrast of structures that can be expected in ACDM with 
gaussian primordial perturbations, nor does it account for the fact that not all structures will 
be identical. By contrast our analysis models the actual expected distribution of void depths 
and sizes. Although, as in our analysis, all structures are assumed to be centred at the same 
redshift for the purposes of calculating the template ISW shift for individual structures, when 
lines of sight overlap their template contributions are simply added. 

In this way Papai et al. obtain the template ISW temperature map, scale this template 
map by a factor A, and maximise the likelihood function 



where Tcmb is the WMAP ILC map, Tisw is the constructed template ISW map and C is 
the pixel covariance matrix obtained from best-fit ACDM. 

The value of A that maximises the likelihood L(X) can be interpreted as the average 
central density contrast in the superstructures [36] 



where superscripts 50c and 50t> refer to the average of 50 superclusters and 50 supervoids 
respectively (note that in the linear approximation voids and clusters are entirely equivalent 
except for a minus sign). This is merely the density contrast that the superstructures must 
have in order for the template to fit the data. 

It is clear from the figures in Ref. [36] that the best-fit values are A > 1 for a wide range 
of radial scales of the structures. Thus, although the initial scenario, where all structures are 
postulated to have exactly \Sq\ = 1 , is already extremely unlikely, in order to fit the observed 
CMB the template ISW map must be further scaled up by a factor A > 1. In this case the 
linear approximation used to calculate the template ISW map must break down (indeed it is 
not clear that the approximation can be valid even when \8q\ = 1). Further note that physical 
voids are restricted to Jo > —1 making the interpretation of A > 1 even more problematic. 
Therefore although Ref. [36] find a reasonable fit for their template when A > 1, this value 
is indicative of a problem with modelling the observations. This qualitative conclusion is 
entirely compatible with the more careful quantitative calculation presented here. 

Although this conclusion is not highlighted in Ref. [36] it is clear that correctly inter- 
preted, both results in fact indicate a significant tension between the ACDM model of the 
ISW effect due to superstructures and the observation first reported in Ref. [31]. 

5 Summary and prospectives 

We have calculated the integrated Sachs- Wolfe effect expected in ACDM from ~ 100 /i _1 Mpc 
size structures, using the density profiles predicted by the linear theory of gaussian perturba- 
tions [37] . We find that the most extreme superstructures in the SDSS volume will produce an 
ISW signal of < 2 /iK. This matches well with ISW maps generated from A-body simulations 
of the ACDM cosmology [38]. 



L(X) = — - (Tcmb — ATisw) C 1 (Tcmb — ATisw) 



(4.2) 




) 



(4.3) 
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Our result is about 3 times larger than earlier calculations [33, 34] which assumed 
compensated top-hat density profiles. Such an assumption, while well-motivated for non- 
linear structures that have formed at small scales [40] , should not apply on the large scales of 
the superstructures considered here. Ref. [35] noted that the ISW signal from such structures 
should therefore be larger and our results confirm this. 

Nevertheless we have demonstrated that the ISW signal claimed to have been detected 
in Ref. [31] is still > 3a larger than the signal expected in ACDM. This tension persists 
even after allowing for likely selection effects. In fact, even the most extreme underdensities 
in the SDSS volume would still produce a signal discrepant by > 3cr with the observed 
signal. Therefore the observed signal cannot be due to a selection effect. We concur with 
Refs. [33, 34] that, if the observed signal is due to the ISW effect, deep superstructures appear 
to be far more numerous than expected in a ACDM cosmology. 

This differs from the conclusion of Ref. [36] who also used similar density profiles. We 
have demonstrated that this difference is not due to the small difference in profiles used, and 
have argued that it is instead because those authors interpreted at face value their template- 
fitting method when applied for density contrasts where the template necessarily breaks down 
for voids (i.e., requires S < — 1!). A more appropriate interpretation is entirely compatible 
with the results presented here. 

An interesting question is whether the expected signal of the most extreme super- 
structures in a ACDM universe is possible to detect in principle. For 50 4° circles of 
the CMB, Ref. [31] reports the statistical uncertainty is 3.1 fiK. Hence the expected sta- 
tistical uncertainty in the average for a sample of N such patches will be approximately 
(3.1 x \JhQ/\fN) [iK, so that for N ~ 3000, the uncertainty on a measurement of (AT) is 
AT no ; se ~ 0.4 fjK. Thus a detection of |(AT)| = 1 //K can be made with roughly 2.5a sig- 
nificance, averaged over an ensemble of 3000 superstructures. From Fig. 2 it is seen that for 
the 3000 largest voids in the SDSS DR6 survey volume, |(AT)| ~ 1 /iK which is of the right 
order of magnitude but somewhat too small for detection at high significance. However the 
SDSS window is not large enough to contain 3000 independent 4° patches on the sky so in 
any case a larger survey would be needed in order to measure a statistically significant signal 
and this would contain more supervoids. This order-of-magnitude estimate indicates that 
even if the reported observation is a statistical anomaly or flawed in some other respect, the 
ISW imprint of superstructures in a ACDM cosmology may be large enough to be detected 
in future surveys. 

As noted earlier, the detection of the ISW imprint of individual superstructures pro- 
vides an important complement to full-sky CMB-galaxy cross-correlation studies. It has the 
potential to provide information about the radii, density contrasts and density profiles of spe- 
cific structures in the dark matter distribution that lie in the extreme tail of the probability 
distribution function. Our calculation demonstrates that the predicted ISW signal from the 
most extreme superstructures is far too small to explain the temperature fluctuations seen in 
Ref. [31], and that therefore the reported observations are in tension with the ACDM model. 
Further observational investigation of this issue is warranted. 

Nevertheless, the standard cosmology successfully fits many other observations, and any 
resolution of the discrepancy discussed here must also satisfy these constraints. A possible 
physical explanation would be that the primordial perturbations are non-gaussian. This 
would influence both the abundance of these extreme regions (e.g., [51, 52]) as well as their 
density profile, thus changing their expected ISW signal. A primordial non-gaussianity dis- 
proportionately affects the tails of the probability distribution function. Since the observation 
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in Ref. [31] specifically probes the ISW effect of structures in these tails, it is expected to 
be more sensitive to non-gaussianity than the standard cross-correlation ISW studies. In- 
terestingly, several cross-correlation studies also find an enhancement of the ISW signal over 
ACDM expectations, e.g. [28, 29, 55], though given the large error bars this discrepancy 
has not been considered too significant. Therefore, primordial non-gaussianity may be able 
to explain this signal while preserving the success of ACDM on other fronts. However, this 
would then undermine the use of the ISW effect as an independent test for A. 

It is worth noting that, qualitatively, a similar enhancement is found for overdensities 
as well as for voids, though we did not quantify this. As pointed in Ref. [53] a primordial 
skewness, parameterised by /nl> would not be able to enhance the abundance of both over- 
and underdense regions simultaneously; however a primordial kurtosis, parameterised by a 
positive #nl> would indeed do so, and be less constrained by the CMB. 

Another possible explanation might lie in a modification of the growth rate of pertur- 
bations as can happen in, e.g., models based on scalar-tensor gravity [54]. Such models may 
also change predictions by altering the behaviour of photons in gravitational potentials. The 
presence of large-scale inhomogeneities can themselves alter the growth rate and this too 
deserves further attention. 
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